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We study the entanglement of a two-qubit one dimensional XYZ Heisenberg chain in 
thermal equilibrium at temperature T. We obtain an analytical expression for the con- 
currence of this system in terms of the parameters of the Hamiltonian and T. We show 
that depending on the relation among the coupling constants it is possible to increase 
the amount of entanglement of the system increasing its anisotropy. We also show nu- 
merically that for all sets of the coupling constants entanglement is a monotonically 
decreasing function of the temperature T, proving that we must have at least an exter- 
nal magnetic field in the z-direction to obtain a behavior where entanglement increases 
with T. 

Keywords: Thermal entanglement, spin chains 



1. Introduction 

Since the beginning of Quantum Mechanics (QM) it was recognized that the su- 
perposition principle imphed novel and counter-intuitive experimental results i^^-^. 
Erwin Schrodinger used the word entanglement to succinctly express one of the 
most striking non-local features of that new theory. However, for a long time since 
the discovery of entanglement, nobody expected that it would have some practical 
application and the majority of the physicists thought that it would be restricted 
to the conceptual discussions of QM. Nevertheless, the last decade of the twentieth 
century revealed incredible practical applications for the entanglement. From new 
quantum algorithms that outperform their classical counterparts to quantum 
communication ^'^'^ we see the usefulness of quantum entanglement. 

Once accepted that entanglement is a resource we can manipulate to do useful 
tasks, it became a necessity to quantify it. It turned out that this quantification is 
not an easy task and up to now it is one of the most challenging open questions 
in Quantum Information Theory S'lO-H. Fortunately, for bipartite systems of two 
levels, i. e., a pair of qubits, there exists an analytical expression to quantify its 
amount of entanglement It is called Entanglement of Formation {Ep) ^. Given 
the density matrix p that describes our pair of qubits, Ep is the average entan- 
glement of the pure states of the decomposition of p, minimized over all possible 
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decompositions: 

EF{p)=mmY,p,E{i',), (1) 

i 

where Y^iPi = 1, < < L and p = Y^^Pi \tpi) {tpi\. Here E{tp) is the von 
Neumann entropy of either of the two qubits ^'^ . For a pair of qubits Wootters et al 
have shown that Ep is a monotonically increasing function of the concurrence 
C, which one can prove to be an entanglement monotone. When C = 1 we have 
maximally entangled states and when C = we do not have entanglement. Since 
C is mathematically simpler to deal with than Ep and that knowing it we can 
automatically get Ep, we concentrate our efforts calculating C to study the amount 
of entanglement between two qubits. The concurrence between them is ^^: 

C = max{Ai - A2 - A3 - A4, 0}, (2) 

where Ai, A2, A3, and A4 are the square roots of the eigenvalues, in decreasing order, 
of the matrix R = pp. Here p is the spin flipped matrix given by 

p = (a-y O ay) p* {ay O ay) . (3) 

The symbol p* means complex conjugation of the matrix p in the standard basis 

{|00),|G1),|10),|11)}. 

Recently, considerable efforts have been devoted in the study of Heisenberg spin 
systems I4,i5,i6,i7,i8,i9 concerning their amount of entanglement for some tempera- 
ture T. These systems are relatively simple and can describe real solid state systems 

which can possibly be used to generate entangled qnbits 21,22,23 nsed in the 
construction of quantum gates 24,25^ building blocks of any quantum computer. 

In this article we deal with ID Heisenberg open chains with no external magnetic 
field and only nearest neighbor interactions. This system is often called XYZ model 
and the Hamiltonian that describes the system is: 

i=l ^ ^ 

where we restrict to chains of only two spins {N = 2), Jx,Jy, Jz are the coupling 
constants, a^, ay, a^ are the Pauli matrices, and h= 1. 

Here we present an analytical formula for the concurrence and a detailed study 
for the general XYZ model. We show that there exist regions for the XY and XYZ 
models where the concurrence increases monotonically as we increase the anisotropy 
of the system. This behavior is a new feature of thermal entanglement in chains with 
no external magnetic field, since Wang ^"^ and Kamta et al only studied systems 
where the concurrence decreases with anisotropy. In these regions the anisotropy 
also increases the critical temperature Tc beyond which the concurrence is zero. 
We review the results obtained by Refs. 16-18 for the XY model {J^ = 0) and by 
Ref. 15 for the XXX model {J^ = Jy = Jz), which are obtained as particular cases 
of our general solution. We also numerically show that contrary to the cases where 
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we have an external magnetic field I5'i6,i7,i8^ there does not exist any set of the 
coupling constants which allows an increase of the concurrence as we increase the 
temperature T of the system. 

2. XYZ Thermal State: An Overview 

In order to study the XYZ model we rewrite the Hamiltonian (4) in the following 
form: 



H = ^alal + ^^alal + ^^y^, (5) 



where A = J^ — Jy and S = Jx+Jy The four eigenvectors of the Hamiltonian are the 
four Bell states (maximally entangled states): iJ = A$± |$^) and H\^^) = 
X^± |vI/±), where |<I>±) = (l/v^) (|00) ± |11)), |*±) = {I/V2) (|01) ± |10)), A$± = 
(J^ ± A) /4, and A^± = (- J, ± S) /4. 

We use the parameter 6 — A/E to measure the anisotropy of the system 
When ^ = and = we have the isotropic XY model and when 5 = ±1 and 
J2 = wc have the Ising model. 

The density matrix which describes a system in equilibrium with a thermal 
reservoir at temperature T (canonical ensemble) is p = exp{—H/kT)/Z, where 
Z = Tr {exp (— iJ/fcT)} is the partition function and k is Boltzmann's constant. 
Hamiltonian (5) gives the following thermal state written in the standard basis: 



1 



/ e-"cosh(/3) -e"" sinh(/3) \ 

e" cosh(7) — e" sinh(7) 

— e" sinh(7) cosh(7) 

\-e-"sinh(/J) e-«cosh(/3) / 



(6) 



where a = J^/{AkT), f3 = A/(4fcT), 7 = S/(4A;T), and Z = 
2 (exp (—a) cosh(/3) + exp (a) cosh(7)). 

The square roots of the four eigenvalues of the matrix R = pp, which enable us 
to obtain the concurrence, are: 

Af = ^(cosh(/3)±sinh(/3)), (7) 

A± =^(cosh(7)±sinh(7)). (8) 

It is not a trivial task to put in decreasing order Xf and Xfj, since we need the 
values of a, (3, and 7 to correctly arrange them. However, the concurrence of p, 
Eq. (2), can be written in the following somewhat simple analytical form: 

rmax{0, Ci},if2a>|/3|-|7|, 

lmax{0,C2},if2a<|/3|-|7|, ^''> 



where 



^ ^ e°sinh(|7|)-e-"cosh(/3) 
^ e" cosh(7) + e-" cosh(/3) ' ^ ' 



I 



I 
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^ _ e-"sinh(|/3|)-e"cosh(7) 
^ e« cosh(7) + e-« cosh(/3) ' ^ ^ 

Before we study in detail Eq. (9) for arbitrary values of Jx,Jy, and Jz we will 
analyze the concurrence for some particular interesting cases. 

3. Ising Model 

In the Ising model Jx = Jy = This implies that /3 = 7 = 0. Substituting in 
Eq. (9) we obtain: 



We can understand why the thermal Ising system is not entangled for any T looking 
at the density matrix (6). When /3 = 7 = it is diagonal in the standard basis im- 
plying no quantum correlations. This result is not surprise because despite of having 
four maximally entangled states as eigenvectors. |<t>*) and l^*^) are degenerated, 
which implies that the Ising thermal state has no entanglement. 

4. XY Model 

When Jj, = we deal with the XY model. It is called isotropic XY model when 
Jx = Jy = J and anisotropic XY model when Jx ^ Jy. We study separately the 



4.1. Isotropic Case 

In the isotropic XY model a = (3 = and 7 = J/(2fcT). Using these values in 
Eq. (9) we get: 



Analyzing Eq. (13) we see that for very low temperatures we have concurrence 

close to 1 and that it decreases monotonically with the temperature until a critical 
temperature which is the solution of sinh(| J|/(2fcTc)) = 1. We also see that 

systems with strong coupling (higher J) have a greater concurrence for a given 
T than those systems with weak coupling and that the concurrence is the same 
whether the system is ferromagnetic (J < 0) or antiferromagnetic (J > 0). 

4.2. Anisotropic Case 

Setting = in Eq. (9) we obtain the following expression for the concurrence of 
the anisotropic XY model: 




(12) 



two cases. 




(13) 




(14) 
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/ sinh(|/3|)-cosh(7)\ 

C = max < 0, — r-r^ > , if d > 1. (15) 

\ cosh(7) + cosh(/3) j ' i i - ^ ' 

Looking at Eqs. (14) and (15) we see two regions of anisotropy. The first region, 
1^1 < 1, was studied by Wang and by Kamta and Starace They have shown 
that increasing the anisotropy parameter \S\ the concurrence decreases for a given 
temperature T and that when = 1 the concurrence is zero for all T (Ising model). 
However, in the second region, \5\ > 1, we see that the concurrence increases as we 
increase the anisotropy parameter 6 and that the critical temperature Tc increases 
as we increase the anisotropy. Wc also sec that in both regions the concurrence is a 
monotonically decreasing function of the temperature. Figs. (1) and (2) show this 
behavior. 




Fig. 1. The dependence of the concurrence C for the XY model with the absolute temperature 
kT. The solid line represents 5 = 0.3, the dotted line shows 5 = 0.6, and the dashed line is for 
S = 0.8. We clearly see that the greater S the lower is C. We have set S = 1. 




Fig. 2. The dependence of the concurrence C for the XY model with the absolute temperature 

kT. The solid line represents <5 = 1.2, the dotted line shows S = 1.4, and the dashed line is for 
(5 = 1.7. We clearly see now that the greater 5 the greater is C. We have set S = 1. 



I 



I 



February 1, 2008 7:19 WSPC/INSTRUCTION FILE IJQIl 



6 G. Rigolin 

It is interesting to note that the conditions on \S\ given in Eqs. (14) and (15) are 
equivalent to J^Jy > and J^Jy <0 respectively. This implies that the anisotropy 
increases the concurrence if and Jy have different signs. 

5. XXX Model 

When = Jy = = J we have the XXX model , which implies that /3 = and 

7 = 2a. Eq. (9) then shows that C = if J < and 

C = max|0, ^-p^^j, if J>0, (16) 

It is interesting to note that contrary to the XY isotropic model the concur- 
rence for the fcrromagnc'tic XXX model is always zero This can be understood 
if we consider the eigenvectors of the system. For J < we have a degener- 
acy for the ground state, which is formed by the triplets. Therefore p{T = 0) 
= (1/3) (|^'+) (^-+1 + |$+) ($+1 + 1$-) ($"|), which is a non-entangled state. And 
increasing the temperature the singlet mixes with the triplets in a non-entangled 
state. However, when J > 0, the ground state is the singlet, a maximally entangled 
state, and we have C = 1. Increasing the temperature reduces the entanglement 
because we mix the triplets with the singlet. See Fig. (3). 

1 

. 8 
0.6 

u 

0.4 
0.2 

0.3 0.6 0.9 1.2 1.5 

kT 

Fig. 3. The dependence of the concurrence C with the absolute temperature kT for three values 
of the coupling constant J for the XXX model. The solid line represents J = 1.5, the dashed line 
shows J = 1, and the dotted line is for J = 0.5. 

6. XXZ Model 

If Jz ^ J and Jx = Jy = J we have the XXZ model. Now /? = and Eq. (9) gives 
C = Oif 2q!< -|7| and 

c = max (o, ^^;^i"MM)- n .^^^^ 

\ e2« cosh(7) -I- 1 J 
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Again we can understand why for 2a < — I7I there exists no entanglement even at 
T = by fooking at the ground state of p. In this region < and |<f>*) are the 
degenerated ground states. Therefore, p{T = 0) = (1/2) (|$+) ($+| + |$") ($"|) 
= (1/2) (|00) (00| + |11) (11|), which is a separable state. By increasing T we mix 
|\E'^) with 1$^), producing a non-entangled state. On the other hand, when 2a > 
— I7I the ground state is the singlet, a maximally entangled state. By increasing 
the temperature the components of the triplet mix with the singlet decreasing the 
concurrence. See Fig. (4). 




Fig. 4. The concurrence C as a function of the absolute temperature kT and of J. We have set 
= —0.5. It is clear that there exists a region where C = for any T. 

It is worthy of mention that we have numerically observed that when > 0, 
an increase of its value increases the entanglement for a given T. We also note 
that whatever the sign of J, increasing its modulo always increases entanglement. 
These two behaviors can be understood analyzing the eigenvalues of the system: 
increasing the absolute value of J or increasing Jz increases the proportion of singlet 
(J positive) or the proportion of I'l'^) (J negative) in the thermal state. These two 
facts are responsible for an increment in the concurrence. 

7. XYZ Thermal State: Detailed Study 

We now analyze the XYZ model. Jx, Jy, and Jz are allowed to assume any values 
and we must study Eq. (9) in its general form. We first observe an interesting 
situation. Whenever 2a = \(3\ — I7I we have null concurrence even for T = 0. This 
condition is equivalent to 2Jz = |A| — |S|. This implies that systems with coupling 
constants near this region are not useful in generating entanglement. Figs. (5, 6) 
and (8, 9) highlight this behavior. 

Eq. (9) implies that for fixed Jz there exist regions where increasing the 
anisotropy parameter 5 we increase the concurrence. In the region where 2a < 
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Fig. 5. The concurrence C as a function of A and of kT. We see that there are regions where an 
increase in the anisotropy increases C and that C is a monotonically decreasing function of kT. 
S = 2 and = 1. 
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Fig. 6. The concurrence C as a function of A for various values of kT. For the solid line kT = 0.05, 
dotted line kT = 0.1, short dashed kT = 0.2, dashed kT = 0.4, and long dashed kT = 0.8. We 
clearly see that for | A| > 4, i. e. {2a < — |7|), the higher the anisotropy (greater |A|) the more 
entangled is the system. S = 2 and Jz = 1. 

\(3\ — I7I the higher the anisotropy the more the system is entangled. However, in 
the region where 2a > \(3\ — I7I, which reduces to the region studied by Refs. 16 
and 17 if we put Jz = 0, the anisotropy decreases entanglement. Figs. (5) and (6) 
illustrate this. 

We can understand the behavior of C as we vary A looking at the probability 
distribution P of the four eigenvectors of the XYZ Hamiltonian in the thermal 
state. Here, 

P,.-Tr{\<,^){<,^\p} = '-^^i^^fII^, (18) 
P,.=Tr{\^^){^^\p} = '-^^i^^fI^. (19) 
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We see that when 2a — — \"f\ the thermal state is an equal mixture of |$^) 
(<i>^ for /3 > and (f>+ for /? < 0) and plus a tiny fraction of |^'+). 

The density matrix that describes the system can be written as p — (1/2 — 
e/2) + I*-) {^-\) + e|*+) (*+|, where e < 1. This density matrix is 

separable for e < 1/2, which explains why in the region where 2a — \P\ — I7I we 
have no entanglement. See Fig. (7). 
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Fig. 7. The probability distribution P of the eigenvectors of the XYZ Hamiltonian in the thermal 
state as a function of A. The solid/red line gives P for dotted/blue for I"!"") , short 

dashed/black for and long dashed/green for |^~)- For large |A| and |A| !5J only one 

maximally entangled state dominates, justifying why we have high concurrences in this region. 
And more, as we increase the temperature, we see that near A = the triplets mix with the 
singlet decreasing C. For large |A| this mixing only occurs for higher temperatures. Jz = 1 and 
S = 2. Top: kT = 0.4. Bottom: kT = 0.8. 

We can also understand why for 2a < \(3\ — \^\ the anisotropy increases en- 
tanglement studying Eq. (18) and (19). In this region, we see that an increase in 
the anisotropy (higher A) produces a thermal state almost dominated by only one 
maximally entangled state, which implies an increase in C. On the other hand, if 
we are in the region where 2a > \f3\ — I7I, an increase in the anisotropy produces 
a mixed state of two maximally entangled state, which causes a decrease in C. 
In the limit where 2a = |/3| — I7I we have a mixture of these two states in equal 
proportions, implying C = 0. See Fig. (7). 
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Fixing Jx, Jy, and varying we see that the concurrence increases if we pick 
values of Jz greater or lower than (|A| — |I]|)/2. There exists, however, a value of 
Jz beyond which C does not increase anymore. This behavior is more drastic if we 
are in the region where 2a > — |7|. There, only for kT « we obtain C « 1. 
For any other value of kT, increasing Jz makes C — > Cmax, where Cmax < 1- And 
more, the higher kT the lower is the value of Cmax ■ If we are in the region where 
2a < \p\ — I7I, decreasing Jz we still can asymptotically obtain C — 1 for kT > 0. 
See Figs. (8) and (9). 




Fig. 8. The dependence of the concurrence C as a function of Jz and kT. As we move away from 
Jz = 3, i. e. (I A| — |S|)/2), we get higher values for C, which is a decreasing function of kT. A = 7 
and S = 1. 
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Fig. 9. The concurrence C is a function of Jz for different values of kT. For the solid line 
kT = 0.05, dotted line kT = 0.1, short dashed kT = 0.2, dashed kT = 0.4, and long dashed 
kT = 0.8. In the region where 2a > — I7I increasing Jz we get C Ri 1 only for kT f» 0. A = 7 
and E = 1. 



We can again get a physical picture of the behavior of C as a function of Jz 
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studying the probability distribution P of the four eigenvectors of the XYZ Hamil- 
tonian in the thermal state. We see that as we move away from — (|A| — |S|)/2 
one of the four maximally entangled states begins to dominate, explaining why C 
increases. But only in the region where 2a < \f3\ — I7I there exists for kT > a 
value of \Jz\ beyond which the probability distribution is zero for the other three 
Bell states, implying C = 1. If we are in the region where 2a > — I7I we get 
for large a reasonable contribution of another maximally entangled state (|^'"*')), 
justifying why C < 1 in this region. See Fig. (10). 
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Fig. 10. The probability distribution P of the eigenvectors of the XYZ Hamiltonian in the thermal 
state as a function of Jz. The solid/red line gives P for l"!"^), dotted/blue for , short 

dashed/black for IvE"^), and long dashed/green for |^~). As we move away from = 3 one of 
the four maximally entangled states begins to dominate, explaining why C increases. A = 7 and 
S = 1. Top: kT = 0.4. Bottom: kT = 0.8. 



We end our study of the XYZ chain pointing that we have numerically searched 
for a set of coupling constants that should give dC/d{kT) > 0. If such set existed 
we would have found a region where increasing the temperature causes an increase 
in the entanglement with no need of an external magnetic field I4,i5,i6^ -yyg tested 
the sign of dC/d{kT) for all combinations of Jx,Jy, and ranging from —2 to 2 
in increments of 0.01 and from —50 to 50 in increments of 0.5. We have not found 
any set of constants where dC/d{kT) > 0. This result suggests that we must have 
an external magnetic field in order to get a positive derivative. 
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8. Conclusion 

In this article, we have studied the thermal entanglement for the XYZ Heisenberg 
chain. We restricted our attention to chains of two qubits with no external magnetic 
field. 

We reviewed the well known results for the XY and XXX models and analyzed 
the XXZ and XYZ models in detail. 

We presented a general analytical expression for the concurrence of the XYZ 
model and then particularized to some simple cases. We have shown that there 
exist regimes in the XY and XYZ models where an increase of the anisotropy 
in the coupling constants causes an increase in the amount of entanglement for 
a given temperature T. In these regions, the critical temperature, beyond which 
the entanglement is zero, also increases with the anisotropy. This result is very 
interesting since up to now only regions where a decrease of the entanglement with 
anisotropy were studied. 

Wc have numerically checked that there docs not exist any combination of the 
coupling constants where the entanglement increases with T. This result suggests 
that we must have an external magnetic field applied to the qubits to prevent the 
entanglement to monotonically decrease with T, i. e., we need external fields to 
obtain a region where the entanglement increases as wc increase T. 

Finally, we want to mention that we are just begimiing to understand the re- 
lation between entanglement and temperature in Heisenberg chains. Here we have 
dealt with only two qubits and nearest neighbors interactions. It would be interest- 
ing to study chains with a large number of qubits as well as chains with long range 
interactions. 
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